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OPERATING  CHARACTERISTICS  FOR  DETECTION 
OF  A FADING  SIGNAL  WITH  K DEPENDENT  FADES 
AND  D-FOLD  DIVERSITY  IN  M ALTERNATIVE  LOCATIONS 


INTRODUCTION 


The  probabilities  of  correct  decision  and  false  alarm  for  a fading 
signal  with  D-fold  diversity  in  M alternative  locations  were  derived 
and  evaluated  numerically  in  reference  1 for  a wide  range  of  values  of 
D,  M and  the  signal-energy-to-noise-density  ratio.  It  was  assumed  in 
that  report  that  all  the  M decision  variables  were  statistically  inde- 
pendent of  each  other.  In  addition,  the  DM  individual  diversity-branch 
nuts  were  assumed  to  be  statistically  independent  of  each  other. 

In  this  report,  we  want  to  generalize  the  situation  to  include  the 
ility  of  having  K signal  samples  (within  each  diversity  branch) 
at  are  completely  statistically  dependent  on  each  other.  A possible 
communication  scheme  that  could  lead  to  this  circumstance  is  illustrated 
in  figure  1.  For  this  example,  there  are  two  alternative  locations 
(M  = 2)  for  the  signal,  and  there  are  three  separated  tones  transmitted 
simultaneously  (D  = 3).  However,  instead  of  building  narrowband  filters 
of  duration  T for  each  diversity  branch  at  the  receiver  to  include  all 
the  received  signal  energy  in  each  diversity  branch,  two  narrowband 
filters  of  duration  T/2  are  used  in  each  diversity  branch,  and  their 
envelopes -squared  are  summed.  Although  the  two  noise  segments  (for 
K = 2)  are  independent  (for  white  noise),  because  they  do  not  overlap 
in  time,  the  amplitudes  of  the  signal  fades  can  be  completely  dependent 
if  T is  significantly  less  than  the  correlation  time  of  the  channel. 

It  is  of  interest  to  know  the  degradation  in  performance  caused  by  this 
suboptimum  processing  technique.  The  results  in  reference  1 correspond 
to  K = 1. 

This  report  is  a generalization  of  reference  1.  To  avoid  duplica- 
tion, that  report  is  referenced  for  the  physical  motivation  and  applica- 
tion, background  information,  and  interpretations.  In  addition,  the 
definitions  of  the  error  probabilities  of  interest,  the  decision-making 
scheme,  and  the  signal  and  noise  models  are  documented  in  reference  1. 

In  this  report,  knowledge  of  the  material  in  the  earlier  report  is 
assumed,  and  reference  is  made  to  the  results,  symbols,  and  notation  in 
that  report. 
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Figure  1.  Time-Frequency  Occupancy  Pattern 


PROBLEM  DEFINITION  AND  PERFORMANCE  CALCULATIONS 


The  signal  diversity  is  D;  this  is  the  number  of  statistically 
independent  signal  components  (occupied  frequency  slots)  per  transmis- 
sion. K is  the  number  of  independent  noise  samples  per  independent 
signal  component.  Hence,  DK  is  the  total  number  of  independent  noise 
components  in  one  of  the  M alternative  signal  locations. 

If  a signal  is  present  in  alternative  location  1 (without  loss  of 
generality),  the  M variables,  on  which  decisions  about  signal  absence 
or  presence  and  identity  must  be  made,*  are  given  by  (see  reference  1, 
equations  (24)  and  (B-l)) 

D K 

‘.-2  2 

d=l  k= 1 

D K 

%•  2 2 

d=l  k= 1 


*The  optimum  processor  is  addressed  in  a later  section.  Here  we 
consider  the  simple  sum  of  envelopes-squared  processor. 


sa  * "dk’ 


(n 


(m) 

dk 


, 2 < m < M 
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Here,  is  the  receiver  output  signal  voltage  complex  envelope  on 
the  d-th  diversity  branch;  notice  that  it  is  not  dependent  on  k.  The  D 
random  variables  {s^}  are  all  statistically  independent  of  each  other. 
The  random  variable  {n (PJ  } is  the  output  noise  complex  envelope  in  the 
k-th  sampling  interval  xn  the  d-th  diversity  branch  of  the  m-th  alter- 
native location.  All  Dktl  random  variables  (n^  } are  statistically 
independent  of  each  other  and  of  the  variables  { s ^ } . 

It  follows  from  the  above  that  the  M decision  variables  {xm}  in 
equation  (1)  are  statistically  independent  of  each  other.  If  the  prob- 
ability density  function  of  xm,  2 £ m < M is  denoted  by  pg.and  that  of 
x by  Pj , then  the  probability  of  false  alarm  (see  reference  1,  page  4) 
is 

PpA  = 1 - [P0(A)]M  , (2) 

and  the  probability  of  correct  decision  is 


P CD  = / dx  (x) [P0(x)]M-1  , (3) 

«/A 

where  A is  a threshold. 

From  reference  1,  equation  (A- 10),  we  have  the  upper  and  lower 
bounds  on  P^: 

M-l 

[1  - PpA]  M [i  . Pi(A)]<  PCD  ll  - P i (A)  . (4) 

Since  P <<  1 for  our  applications,  this  is  an  exceedingly  useful 
and  tight  bound  for  all  M;  it  requires  only  that  we  be  able  to  evaluate 
the  cumulative  probability  distribution  Pi,  rather  than  perform  the 
numerical  integration  indicated  by  equation  (3). 

The  derivation  of  Pi  is  given  in  appendix  A.  Two  alternative 
forms  are  given  by 


I - PiCA)  - exp(-A)(l  - r)D5;(D  + ; - 

n=0  ' 

= l-CI-V)DJ//D  + nn-AYn[1.exP(-A)eDK_1+n(A)] 

n=o'  ' 


(5) 
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Reasons  for  choosing  one  form  over  the  other  and  a method  for  accurate 
evaluation  of  equation  (5)  for  large  D are  given  in  appendix  A.  The 
parameter 


where 


R = 77  =- 


is  the  signal-to-noise  power  ratio  of  an  individual  sample  in  equation 
(1) . The  quantity 


? (A)  = ^Ak/k!  = 1 + A + ...  + A^/j  ! 


is  the  first  j+1  terms  of  the  power  series  expansion  of  exp(A). 


For  zero  signal-to-noise  ratio,  R = 0,  we  obtain  from  equation  (5), 


1 - P0(A)  = expf-A)eDK1 (A) 


and  then  equation  (2)  yields  the  false  alarm  probability. 


PFA  = 1 - [1  - expC-A)eDK1(A)  ] 


Given  values  of  DK  and  M and  a desired  value  for  PpA>  equation  (10)  is 
solved  for  A.  This  value  is  substituted  in  equation  (5),  which  is  then 
evaluated  as  a function  of  R. 


The  parameter  R in  equation  (7)  can  be  Interpreted  differently  for 
the  model  where  the  received  signal  is  narrowband  and  deterministic 
(except  for  phase)  and  subject  to  slow  Rayleigh  fading.  This  considera- 
tion is  taken  up  in  appendix  B,  with  the  result  that 
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and 


VN° 

~KD 


(11) 


El  VN0 
Nq"  = ~D 


(12) 


where 

Es  is  the  average  received  signal  energy  on  one  diversity  branch 
during  one  individual  processing  interval  (D  = 1,  K = 1), 

E-p  is  the  total  average  received  signal  energy  over  all  diversity 
branches  and  processing  intervals, 

Ej  is  the  average  received  signal  energy  on  one  diversity  branch 
during  all  K processing  intervals,  and 

N0  is  the  (single-sided)  noise  density  level. 

The  curves  to  follow  will  employ  the  signal-energy-to-noise-density 
ratio  parameters  in  equations  (11)  and  (12). 


PERFORMANCE  RESULTS 


In  figures  2_  through  6,  the  probability  of  correct  decision  P„  is 
plotted  versus  Ej/Nq  in  decibels,*  with  diversity  D as  a parameter 
(solid  curves)  for  M = 1,  K = 2,  and  Pp.  = 10_n,  n = 2,  3,  4,  6,  8, 
respectively.  The  dashed  curves  in  each  figure  connect  points^  of  equal 
average  total  received  signal-energy-to-noise-density  ratio,  E^/Ng. 

In  addition  to  the  relevant  observations  made  in  reference  1,  page 
9,  the  following  points  are  worth  noting.  For  Pp.  = 10~2  in  figure  2, 
the  difference  in  minimum  E_/Ng  required  versus  tne  K = 1 case 
(reference  1,  figure  1)  is  1.1  dB  over  a wide  range  of  Pc„  (above  0.9). 
Also,  for  a specified  value  of  Prn,  the  optimum  order  of  diversity  is 
the  same  for  K = 2 a s for  K = 1.  ^However,  both  Ej/Nq  and  E /Nq  must 
be  1.1  dB  greater.  T 


♦That  is,  10  log  (Ej/Ng). 
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For  PpA  = 10'8  in  figure  6,  the  difference  in  required  Eif/No  between 
the  K = 2 and  K = 1 cases  is  1.0  dB  over  a wide  range  of  PcD-  Thus,  the 
difference  in  required  signal -to-noise  ratio  is  relatively  independent 
of  PpA>  at  least  for  M = 1. 

In  figures  7 through  11,  the  curves  for  M = 16,  K = 2 are  plotted. 

We  find  that  the  difference  in  required  signal-to-noise  ratio  versus  the 
K = 1 results  (reference  1,  figures  16  through  20)  is  approximately  1.0 
dB  over  a wide  range  of  P^p  and  Pp^. 

Continuing  on  to  M = 256,  K = 2 in  figures  12  through  16,  we  find 
a difference  in  required  signal-to  noise  ratio  of  approximately  1.0  dB 
(0.9  dB  for  PpA  = 10’8)  between  these  results  and  those  for  K = 1.  Thus, 
the  degradation  in  performance  is  relatively  insensitive  to  the  values 
of  M,  P^q,  and  PpA. 

The  results  for  K = 3 are  presented  in  figures  17  through  31,  again 
for  M = 1,  16,  and  256.  Once  again,  the  same  conclusions  on  invariance 
result,  except  that  the  degradation  now  is  larger,  1.7  dB  for  M = 16, 
for  example. 

Finally,  for  M fixed  at  16,  curves  for  K = 4,  5,  6,  8,  and  10  are 
presented  in  figures  32  through  56.  They  enable  us  to  extract  the 
degradation  behavior  with  change  in  K that  is  plotted  in  figure  57. 
Namely,  for  M = 16,  we  find  that  additional  signal-to-noise  ratio,  to 
the  extent  of  approximately 


3.5  log (K)decibels,  (13) 

is  required,  relative  to  the  K = 1 case,  to  maintain  the  same  perfor- 
mance, PpA,  P(-D.  This  is  a rough  rule-of-thumb;  the  exact  values  can 
be  deduced  from  the  curves  in  figures  2 through  56  or  by  modifying  the 
program  in  appendix  A to  fit  the  user's  case. 
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OPTIMUM  PROCESSOR 


.With  dependent  fading,  the  optimum  processing  of  random  variables 


(x*:  (where  we  had,  for  the  case  in  equation  (1) 


rO)  _ 


I (l)  | 2 

Sd  + "dk  ' m 


. Cm)  1 2 


, 2 < m < M 


is  not  simply  to  sum  them.  Rather,  we  find  in  appendix  C that,  for 
K = 2,  we  should  consider  the  test 


D / v D 

maX  V (A™)  * rt("0\  * V 0nr  1 

< m < M Zr d!  + ad2  ) + Z*nM 
d=l'  ' d=l 


v di  d2  , 


where 


. Cm)  _ 


1 + 2R 


can  be  interpreted  as  a scaled  and  normalized  envelope -squared  sample. 
Although  equation  (15)  is  not  much  more  difficult  to  implement  than  the 
usual  sum  (which  is  proportional  to  the  leading  sum  in  equation  (15)), 
the  test  (equation  (15))  requires  knowledge  of  the  absolute  levels  of 
a|  and  ofi.  Also,  the  exact  performance  characteristics  of  equation  (15) 
would  be  impossible  to  determine  analytically. 


' (m)  (m)  , 

a)  « 1, 

di  d2 


IQ(x)  = 1 + %x2,  tnIQ(x)  = \x2  , for  small  x, 
the  second  sum  of  equation  (15)  is  approximately 


r 
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km?  « i(«ff  * “i?)  • t>« 

a=i  d=ix 

Thus,  in  case  equation  (17)  is  true,  the  processor  (15)  takes  the  form 
of  a sum  of  two  squared  envelopes,  which  are  then  summed  over  the  diver- 
sity branches. 

For 


(m)  (m) 
di  ad2 


» 1 


(20) 


since 


Io  Cx) 


exp(x) 


V2ttx 

the  processor  (15)  takes  the  form 

D 


, £nl0(x)  * x,  for  large  x. 


(21) 


d=l ' ' 


(22) 


which  is  a sum  of  envelopes  that  are  then  squared  and  summed  over  the 
diversity  branches.  The  likelihood  of  equation  (17)  or  (20)  occurring 
is  discussed  in  appendix  C.  It  is  shown  that  if  R « 1,  the  processor 
rule  on  the  right  side  of  equation  (19)  is  substantially  optimum, 
whereas  if  R >>  1,  the  processor  rule  in  equation  (22)  is  substantially 
optimum. 


PARTIALLY  DEPENDENT  SIGNAL  FADING 


When  the  rate  of  signal  fading  is  neither  faster  nor  slower  than 
the  rate  of  signaling  (time  T in  figure  1),  but  the  two  are  comparable, 
the  signal  complex  amplitude  s^  in  equation  (1)  cannot  be  considered 
constant  over  all  K samples.  If  we  let  the  signal  amplitude  obey  Ray- 
leigh fading,  we  find  (appendix  D)  that  the  characteristic  function  of 
the  decision  variable  with  signal  present  is  given  by 


K 

nu 

k=  1 


U(on  + 


2x[x))} 


-D 


(23) 


8 


and  where  is  the  real  part  of  the  Rayleigh-fading  signal  complex 
amplitude  s<j  at  the  k-th  sampling  instant. 

For  completely  dependent  fading,  the  elements  of  the  matrix  in 
equation  (24)  are  all  equal  to  a2  = o2/2,  and 


Ko2/2 , k = 1 
0,  1 < k 


(dependent  fading) , 


(25) 


in  which  case  equation  (23)  reduces  to  equation  (A-13).  On  the  other 

hand,  for  completely  independent  fading,  the  matrix  in  equation  (24)  is 

^o2I,  and 
s 

X,W  = *so2,  all  k (independent  fading),  (26) 

1C  s 

DIC 

in  which  case  equation  (23)  reduces  to  [1  - i£;(o^  + °i)]  • This 

corresponds  to  the  case  studied  in  reference  1 for  K = 1 (see,  for 
example,  reference  1,  equation  (B-3)), 


For  signal  absent,  the  characteristic  function  in  equation  (23) 
reduces  to  the  same  result  as  that  to  which  equation  (A-13)  would 
reduce.  Therefore,  the  false  alarm  probability  in  equation  (10)  is 
applicable  to  this  case  as  well. 


A closed  form  for  the  probability  density  function  corresponding  to 
the  general  result,  equation  (23),  is  possible,  but  would  be  extremely 
tedious.  Perhaps  the  best  numerical  procedure  is  straightforward 
Fourier  transformation  of  equation  (23)  (see  reference  2,  for  example). 
But,  in  any  event,  the  results  of  this  report  and  reference  1 furnish 
bounds  on  performance. 
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SUMMARY 


The  numerical  results  here  for  completely  dependent  fading,  coupled 
with  the  results  for  independent  fading  in  reference  1,  serve  to  bound 
the  performance  in  practical  cases  of  partially  dependent  fading.  This 
would  be  true  whether  the  dependence  is  intentional,  as,  for  example, 
sampling  frequently  in  time  in  figure  1 (large  K)  or  unintentional,  as, 
for  example,  not  being  aware  of  the  true  rate  of  signal  fading.  If  the 
losses  indicated  by  figure  57  are  too  large  for  some  larger  values  of  K, 
it  may  be  necessary  to  resort  to  the  more  exact  result  in  equation  (23), 
if  the  fading  is  fairly  independent,  to  obtain  a more  exact  result.  The 
gains  expected  by  doing  optimum  combination  of  individual  samples  is 
probably  minimal  and  requires  knowledge  of  signal  and  noise  levels  that 
often  is  not  available. 
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Figure  4.  Detection  Characteristics  for  K = 2,  M = 1,  P 
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Figure  6.  Detection  Characteristics  for  K = 2,  M = 1,  P_,  = 10 
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Figure  14.  Detection  Characteristics  for  K = 2,  M = 256,  P = 10 
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Figure  20.  Detection  Characteristics  for  K = 3,  M = 1,  P ^ = 10“6 
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Figure  22.  Detection  Characteristics  for  K = 3,  M = 16,  P = 10 
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Figure  23.  Detection  Characteristics  for  K = 3,  M = 16,  P 
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Figure  25.  Detection  Characteristics  for  K = 3,  M = 16,  PpA  = 10 
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Figure  28.  Detection  Characteristics  for  K = 3,  M = 256,  P = 10 
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Figure  31.  Detection  Characteristics  for  K 
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Figure  32.  Detection  Characteristics  for  K = 4,  M = 16,  P = 10 
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Figure  34.  Detection  Characteristics  for  K = 4,  M = 16,  PpA 


10 


-4 


4 3 


TR  5739 


47 


TR  5739 


Figure  41.  Detection  Characteristics  for  K = 5,  M = 16,  t>  = 10 
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Figure  43.  Detection  Characteristics  for  K 
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Figure  44.  Detection  Characteristics  for  K = 6,  M = 16,  P 
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Figure  45.  Detection  Characteristics  for  K = 6,  M = 16,  P = 10~6 
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Ej/N0  (dB) 

Figure  50.  Detection  Characteristics  for  K = 8,  M = 16,  P = 10~6 

FA 


59 


TR  5739 


Figure  51.  Detection  Characteristics  for  K = 8,  M = 16,  Pp^  = 10 
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Ej/N0  (dB) 

Figure  56.  Detection  Characteristics  for  K = 10,  M = 16,  P = 10*8 
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Appendix  A 


DERIVATION  OF  P:(A) 


From  equation  (1) , we  have 


= ^ |5d+  °<d  . 


d=i  k*i 


where 


°Cj  - 


i «'  l2 

I*4"* l 


The  D random  variables  {a <j}  are  statistical ly  independent  of  each 
other.  Our  approach  will  be  to  find,  first,  the  characteristic  function 
of  aj,  then  the  characteristic  function  of  x[ , and,  finally,  the  proba- 
bility density  and  cumulative  distribution  of  Xj. 

We  start  by  finding  the  characteristic  function  of  conditioned 
on  a fixed  known  value  of  s^. 


C/'fU)  = exj> 

= expfif  ^ + I7) 

' nY 

= evp(t  f|Si  + ndkT) 


using  the  independence  of  (n^)}  and  their  identical  statistics, 
n^j^  is  assumed  complex  Gaussian  with  independent  real  and  imagi 


tics.  Now 
imaginary 


parts : 
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nJK  - X + iy, 

X - y = o,  *3  = yl  z 0 

h!p  * 7 + 7 = 2<rj  h 


Then  the  bracketed  term  in  equation  (A-3)  is  given  by 


I'd*  (2ttV)  expj-  *^r  e xf  ('  ^(5dr  + x)  + 1 <?(5<),+  ^ ) 


^i-  i 2o -’t)  exf  ^ 


>fh,r 

I- 1 e2<r’ 


where  we  have  used  the  result 


jdx  evf  (-«  x%/3v)  --  (j)  exj>(-^)  , fie  ex'  > O (A- 

For  later  use,  the  probability  density  function  corresponding  to 
equation  (A-5)  is 


''r(-  J¥i^iL)I°('  ) 


, x > 0, 


which  is  recognized  as  the  probability  density  function  corresponding 
to  the  square  of  a Rician  variate;  recall  that  this  holds  for  K = 1. 

The  characteristic  function  in  equation  (A-3)  is  now  available  by 
using  equations  (A-5)  and  (A-4) : 


fjTK)  = «r(TTT^) 


0 f the  signal  component  sj  in  equation  ( A- 2)  were  different  on  each  of 
the  K processing  intervals,  equation  (A-8)  would  be  replaced  by 
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(A-9) 


We  will  find  use  for  this  result  later.) 


To  find  the  characteristic  function  faj(£),  we  have  to  average 

equation  (A-8)  over  the  statistics  of  s^.  We  assume  that  sj  is  complex 
Gaussian  with  independent  real  and  imaginary  parts;  then  | s^ | is  a 
Rayleigh  variate.  Letting,  as  in  equation  (A-5) , 


S,  = S,  + . Sd|  , 


$ay  ~ Sdi  = 

- o.  • < A. 


(A- 10) 


is„r  = < , 


we  have  the  average 

= Ifd^  dy  (rr-r;)' evp(~ 

CA- 11") 


If  v is  identified  as  i C K/ (1  - i Ccrjh,  the  average  on  sd  in  equation 
(A-8)  is  now  available  as 


y?)  - 1 fa) 

a d 


[l-'Tu;]  [i-iffy  iKffs1)] 


CA- 12) 
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e characteristic  function  of  xi  is  immediately  available  in 
equation  (A-l),  using  the  statistical  independence  of  {aj},  as 


■ 


(A-l  3") 


r -iD(H  ')  r 

K*V)] 


Since  we  know  the  probability  density  function-characteristic  iunc- 
tion  pair, 

1 , » X1  ,»  0. 

(j-ifo)’  b-')!  “ 


(A- 14) 


we  can  write  the  probability  density  function  corresponding  to  multipli- 
cative characteristic  function  (A-13)  as  the  convolution 


p.w*I 


UJ  gypf-u/a)  (x-u)  expf~  Xb 


du 


U-0!  a 


(p-ft!  b* 


. x > 0, 


(A- 15) 


where 


J s dIk-i),  a = ^ ^ + K 


(A- 16) 


Now,  by  employing  reference  3,  3.383  1 and  8.384.1,  equation  (A-15) 
becomes 


D*  J- 1 


r n exp(->/a)  - 

P'M"  .V(»oTSi 


(D;T  + D;  >(i  - V)),  v>Q.  (A- 1?) 


If  we  scale  random  variable  x according  to 


i _ * « i 

" a <rj 


(A-l 8) 


(we  can  absorb  the  scaling  in  threshold  A)  we  get,  for  the  probability 
densitv  function  of  t. 


A -4 
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* 

?U\=  e*H~- -)  \ 1Pi('D;DK,Yt) 

f’lt)  (DK-O*  (M-KK^  ”V 


exI  Lj  + HR r f>0, 

\ i ( \t>  i'i  V ' ’ 


(A- 19) 


(dk-i)!  (if  kt?) 


where  we  have  used  equation  (A-16)  and  reference  3,  9.212  1,  and  defined 


\ - £ • * = 


KK 
+ KK 


(A - 20 ) 


As  a check  on  equation  (A-19),  we  have  the  following:  tor  K 1, 
equation  (A-19)  reduces  to  equation  (B-3)  in  reference  1;  and  for  R = 0, 
equation  (A-19)  reduces  to  equation  (B-2)  in  reference  1, 

To  find  the  cumulative  distribution  function  of  variate  t,  we 
integrate  equation  (A-19) : 


= J df  f.li') 

Jv 


——2 — [ d*  e.r(-q  tDK  ',f.  (» • » .n). 

(PK-O  ! 0 + KK)  JL 


(A-21) 


By  employing  the  power  series  expansion  of  j F j in  reference  3,  9.210  1 
and  defining  the  partial  exponential 


eiW5  S ^ 


•A"'  +7^ 


(A-22) 


we  find 


The  partial  exponential  in  equation  (A-23)  tends  to  the  constant 
value  exp(A)  as  n -*•  ® and  offers  no  help  in  convergence.  Also,  the 
ratio  of  the  n + 1-th  to  n-th  remaining  factors  tends  to  y as  n -*•  ®. 
Therefore,  equation  (A-23)  is  always  a convergent  series  of  positive 
terms,  since  Y < 1;  see  equation  (A-20) . However,  for  large  KR,  the 
convergence  is  so  slow  that  equation  (A-23)  is  not  very  useful  for 
numerical  evaluation. 

By  taking  the  exp(-A)  term  in  equation  (A-23)  inside  the  sum, and 
adding  and  subtracting  1 from  the  product  of  the  exponential  and  the 
partial  exponential,  we  obtain 


(A-24) 


« = 0 


where  we  used  the  result  (reference  3,  9.100  and  9.121  1), 


2 
h o 


fry  - 


fA-25) 


Now,  the  bracketed  term  in  equation  (A-24)  tends  to  zero  as  n + ®, 
as  do  the  leading  factors.  Thus,  equation  (A-24)  is  a more  rapidly  con- 
vergent series  of  positive  terms  than  equation  (A-23). 

A problem  arises  in  the  accurate  evaluation  of  the  bracketed  term 
in  equation  (A-24) : 


= 1-  exp(-A)  e^^n[A) 

- | “ ^ X r (_  -A  /*"! 

M = 0 

, *„r(-x)  it  A'/"!  • 

' W.  DM* 


(4-26) 
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We  note  the  downward  recursion  of  positive  quantities 

Qn  Z On.,  + C«  , OS 


(A-27) 


where 


C„s  exp(-^APK+y(DKi-n)!  , Osn, 


CA-28) 


and  for  which 


C * exr(-A)-VI*/(CK)!  i C.  = C„,  A/fDK  + n),  |.n. 


iN-1 


Thus,  if  we  evaluate  (Cn}‘  and  Qjj,  we  can  then  evaluate  {Qn) 
via  the  downward  recursion  in  equation  (A-27). 


(A-29) 

N-l 


N 


In  order  to  accomplish  this,  we  evaluate  { C } 1 by  the  forward 

0 


recursion  in  equation  (A-29),  and  then  note  that,  from  equations  (A-26) 
and  (A-28), 


Qkj  = ^K  + KJ^|  ~ (dk+w+i)- 

Jl 


= C 


(A-30) 


n w = 0 (DK  + N + i\ 


This  last  sum  is  terminated  for  desired  accuracy  in  Q^. 

As  a check  on  accuracy,  we  note  from  equation  (A-26)  that 
PK-i 


Q 4-  2exp(-jv)A  /« 1 - I. 


I»^0 


(A-31) 


But  the  terms  in  this  sum  are  computed  in  the  initial  evaluation  of  C0 ; 
so,  if  these  terms  are  saved  and  combined  with  the  final  result  of  Qg 
obtained  via  downward  recursion  (A-27),  the  sum  of  1 dictated  by  equation 
(A-31)  should  be  obtained. 
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The  procedure  is  summarized  below: 

1.  C0 

2.  q,  C2,  ....  CN 


5-  ^N-i’  ^N-2’  •••’  Qo 


At  all  times  we  are  dealing  with  positive  numbers,  and  no  differences 
are  formed.  The  only  difference  occurs  in  the  final  subtraction  from  1 
in  equation  (A-24). 

A FORTRAN  program  for  the  evaluation  and  plotting  of  equation  ( A-24 1 
by  means  of  the  procedure  above  is  presented  on  pages  A-9  through  A- 14. 
The  accuracy  of  the  check  formula  (A- 31)  is  typically  16  decimals  for 
the  double  precision  routine  used. 
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CALL  LINLbG(Z» 1 *20. . YNORM( J) ) 
CALL  SETSMG(Z*30*2. ) 

CALL  L I NtSG ( Z * 9 * 0 . » YNORM ( 1 ) ) 
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Appendix  B 


ALTERNATIVE  INTERPRETATION  OF  R 


Let  the  transmitted  waveform  on  one  diversity  branch  during  one 
individual  processing  interval  (D  = 1,  K = 1)  be 

T>(4)  cosjo^i  + 6 (f\J  = Re  ^ E ft)  e \ . ^B' 1 ^ 

Then  the  received  waveform  corresponding  to  this  transmission  is 

*e{[ftE(O  + N(01  fB'2^ 


where 


= A exj»(i  4>)  (B-3) 

is  the  complex  fading  of  the  channel  (assumed  slow  and  frequency  nonse- 
lective)  and  N(t)  is  the  complex  envelope  of  the  received  noise.  The 
noise  correlation  is 


N (4,)  N*(V)  - 2 N0  £ (4,-  O,  (B-4) 


which  corresponds  to  white  noise  of  single-sided  density  level  N0. 

The  received  signal  energy  in  one  diversity  branch  during  one 
processing  interval  is 


Es  sfdt{AB(t)fs[H+  + elO+  <fl]2  = 4-1^  BJR), 


and  the  corresponding  average  received  signal  energy  is 


V 


OB-5) 


(B-6) 


The  received  signal  energy  on  one  diversity  branch  during  all  K 
processing  intervals  is 


TR  5739 


since  the  fading  is  completely  dependent  during  this  time.  Hence,  the 
average  received  signal  energy  on  one  diversity  branch  during  all  K 
processing  intervals  is 


(B-8) 


The  total  received  signal  energy  over  all  diversity  branches  and 
processing  intervals  is 

£t  = E_,  + E,+  + , CB-9) 

with  average  value 

B,  - F,T>  * Fs  KD.  (B- 10) 

Tins  is  the  average  total  received  signal  energy. 

For  a synchronized  matched  filter  to  the  transmitted  signal  on  one 
diversity  branch  and  one  processing  interval  (D  = 1,  K = 1),  the  output 
complex  envelope  at  the  best  sampling  instant  is  (see  equation  (B-2)) 

y = J cH  [& E [+)  + N Ri]  E* (+'1-  (B-u) 

Hence,  the  squared  envelope  output  is  proportional  to 

jvR  + jdf  NfOE*(+)|  =|s4*f  (B- 1 2) 

Therefore  (see  equation  (A-10)), 

< ■ i?  - ]s? 

and  (see  equation  (A-4)) 

= 2N#  jdf|EfOr  = 2IH,  {dt  (B- i 4 ) 

where  we  have  used  equations  (B-3),  (B-l),  and  (B-4) . As  a result, 


- AJ 


(B-13) 


B-2 
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(B-1S) 


using  equations  (A- 20)  and  (B-6) . Combining  this  result  with  equation 
CB- 10) , we  have 


- KK  IX 


(B-16) 


and,  finally. 


Et  K 


= 'R  K. 


( B— 1 7) 


Equation  (B-16)  is  the  average  total  received  signal-energy-to- 
noise  density  ratio,  whereas  equation  (B-17)  is  the  average  received 
signal-energy  per  diversity  branch-to-noise  density  ratio. 


B-3/B-4 
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Appendix  C 


DERIVATION  OF  OPTIMUM  PROCESSOR 


The  a posteriori  probability  of  alternative  2,  given  (squared-enve- 
lope) observations 


(*o 


< d 5 P,  Uki  K, 


hi 


M, 


(C-l) 


is 


Tr  UIR’l)  - ' ° SiS 


(C— 2) 


by  Bayes  theorem.  Using  the  statistical  independence  of  all  the  variables 
in  equation  (C-l),  we  have,  for  l = 0 (no  signal). 


- IT  Ijr.K’)], 


(C-3) 


(C-4) 


where  the  product  is  over  the  range  indicated  in  equation  (C-l) 

On  the  other  hand,  for  i.  >_  1,  we  have 

'inp(«r  ,-.>«)} 

w.*i 

To  make  an  optimum  selection  among  the  M + 1 alternatives,  we  only 
need  compare  the  numerators  of  the  a posteriori  probabilities  in  equa- 
tion (C-2).  They  are  available  from  equations  (C-3)  and  (C-4).  By 
dividing  these  numerators  by  equation  (C-3),  the  comparison  is  seen  to 
be  between  the  quantities 


J ( f.W, } p-WkIJ 


Or  IS-f  5 M. 


(C-5) 


If  the  a priori  probabilities  are  all  equal,  Pr(l) 
Pi'(M),  the  comparison  is  between 


Pr(2)  = 


C-l 
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W rr)  f'Ki  ■■■■'  x«)  _ 

mo  j uk*,’)  ■■■(’.(*£) 


\ i-  A - M 


Ability  to  evaluate  this  expression  depends  on  the  form  of  the  density 

p fx(*)  XW) 

pllxdl  ’ dK 

In  the  rest  of  this  appendix,  we  restrict  consideration  to  K = 2. 
From  equations  (A-7)  and  (A-4) , for  one  sample  of  a squared-envelope, 
we  have 


^ i. 


But  the  probability  density  function  of  r is  available  from  equations 
(A- 10)  and  (A-ll)  as 


~-~p  r - r-> 


Therefore,  using  a shorthand  notation, 

\ (C-9) 

Substituting  equations  (C-7)  and  (C-8)  into  equation  (C-9)  and  employing 
reference  3,  6.633  2,  we  obtain 


1 -U (~1 v (C*10^ 

v 1 + 21?  C/M'+2R  C )’  \v>°- 


Therefore,  the  desired  ratio  is 


I+2R  OV 


1?  X,+  Xa  W ' K 2 VV,V: 


r I.  m 


■2K  «V 


M (C-ll) 


Then  the  right-hand  term  of  equation  (C-6)  is  given  by 


and  its  natural  logarithm  is 


(*d?-*-x£)  + CC-13) 


Finally,  the  comparison  dictated  by  equation  (C-6)  can  be  summarized 


as 


D 


boax  < 

I cT-  \ 


W'  * «£)  + 3;  (2^vT)j  j a , 


(C- 14) 


where 


10 


») 


Jtc 


11-21? 


liL 


CC-15) 


If  the  maximum  on  the  left  side  exceeds  threshold  A,  we  state  that 
signal  alternative  is  present;  otherwise,  we  state  no  signal  present. 
The  value  of  A is  adjusted  to  realize  a specified  value  of  false  alarm 
probability. 


It  is  of  interest  to  know  what  values  the  argument 


dl  ~d2 

takes  on,  because  we  may  be  able  to  simplify  the  Hn  I0(  ) operation  if 
the  argument  is  often  very  small  or  very  large.  To  this  aim,  we  compute 

1 fin  rrr 

the  average  value  of  Xrf2  for  an  alternative  with  signal  present. 

Since,  from  equation  (C-10), 


(C- 16) 


where 


C-3 
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A - 
B = 

C = 


1 +2R 

<* 

1 + * 

1 

I + 2R 

<r« 

(C- 17) 

R 

2. 

I + 2R 

^ * 

we  have 


oo 

yn=  x,  vj  = ex  fW^^))i;(ora).  (c-j 

0 

By  letting  x„  = t2,  equation  (C-18)  is  converted  into 


6*0 

= 2^dx,  d+y^  f A exp^BCx.+  t1))  I^fCV^O 


(C- 19) 


Employment  of  reference  3,  6.631  1,  enables  the  evaluation  of  the 
integral  on  t and  yields 


'T* 

D o u 

Then,  we  employ  reference  3,  7.621  4,  to  obtain 


(C-20) 


(C-21) 


The  last  step  is  by  use  of  equation  fC-17).  Then  reference  3, 
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9.131  1,  yields 


That  is. 


-0^  if  Y\o  Signal  in  1 Ve  J 


(C-22) 


* 


fC(i+ii)F('^-t-ii(i?r)')  ’f  ,;3ml  j" 

at-fcev'n'Crive  X 


(C - 23) 


There  follows  immediately,  by  means  of  equation  (C-15),  that  the 
average  value  of  the  argument  of  £n  IQ  in  equation  (C-14)  is  given  by 

^ if  o signal  m offevnatM/e  X 


1 1+2-R 


z{ 


^ - 
o(6\  °M2 


a l^e-vnaf  i ve  X 


CC -24) 


From  this  result,  we  conclude  that  if  R <<  1,  all  the  mean  values 
in  equation  (C-24)  are  much  less  than  1,  and  the  expansion  in  equation 
(18)  is  valid.  Alternatively,  if  R >>  1,  since  the  F function  approaches 
4/it  as  R + ® (reference  3,  9.122  1),  we  conclude  that  the  mean  value  for 
the  signal-present  alternative  is  much  larger  than  1,  whereas  the  mean 
value  for  the  other  alternatives  is  u/4.  Although  this  latter  value  is 
not  larger  than  1,  we  don't  know  a priori  which  alternative  contains  the 
signal,  and  we  must,  therefore,  go  along  with  the  signal-bearing  alter- 
native and  use  the  approximation  in  equation  (21)  for  all  alternatives, 
for  large  R.  Considering  the  fact  that  for  large  R the  variable  domi- 
nating the  maximum  comparison  in  equation  (C-14)  virtually  always  will 
be  the  signal-bearing  alternative  anyway  (if  signal  present),  the 
approximate  rule  of  equation  (22)  is  virtually  optimum,  for  large  R. 


V ' 
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Appendix  D 


DERIVATION  OF  CHARACTERISTIC  FUNCTION 
FOR  PARTIALLY  DEPENDENT  SIGNAL  FADING 


The  starting  point  for  this  derivation  is  the  conditional  charac- 
teristic function  (A-9)  for  different  signal  samples  on  the  K intervals: 


We  express  collection  {s^}  in  terms  of  their  real  and  imaginary 
parts  according  to 

(D-2) 

k - t tc=» 

and  assume  that  the  real  and  imaginary  parts,  {x^}  and  (yk),  are  samples 
of  two  independent  zero-mean  Gaussian  processes  x(t)  and  y(t).  Because 
of  the  independence,  the  average  of  the  exponential  in  (D  - 1)  can  be 
expressed  as 


r — 12 

= ■(  dx,  d*K  eyf  (l)^XK')  ? (V<>- •■>  K) 

— * > 


where 
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Now,  for  Gaussian  (x^). 


= (?t)  fcjet  cj'  «xf[-iX  C X], 


(D-5) 


where  column  matrix 


and  square  matrix 


X = [*.- 

r nK 

C = [\  \\ 


(D-7) 


Then,  equation  (D-3)  becomes 


i^7  ^0  = 


(D-8) 


Upon  employment  of  the  integral 


^X«r[-1XTMX]  -fc-f/fc  m)'a, 


(D-9) 


(which  follows  from  equation  (D-5)  when  we  remember  that  p must  have  unit 
volume  for  any  positive  definite  matrix  C) , equation  (D-8)  becomes 


W \ . ^ | 

= [<U(l-2pC)J  . 


<M(C-\r)  * TT 


(D- 10) 


CD-11) 
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where  {X^  } are  the  eigenvalues  of  C.  Therefore, 


y ^ 

^(i-2dC)  = (-2p)K  <U(C-£x)  =TT  (D-12) 


Combining  equations  (D-l) , (D-4),  and  (D-10),  we  get 


--  tUh?  (<+<)] 


(D-13) 


Then,  from  equation  (A-l), 


CD-14) 


which  is  the  characteristic  function  of  the  decision  variable  with 
signal  present  in  it.  This  result  includes  equation  (A-13)  as  a special 
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